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Abstract 

We consider a branching system consisting of particles moving according to a Markov family in K'* and 
undergoing subcritical branching with a constant rate V > 0. New particles immigrate to the system ac- 
cording to homogeneous space-time Poisson random field. The process of the fiuctuations of the rescaled 
occupation time is studied with very mild assumptions on the Markov family. In this general setting a 
functional central limit theorem is proved. The subcriticality of the branching law is crucial for the limit 
behaviour and in a sense overwhelms the properties of the particles' motion. It is for this reason that the 
limit is the same for all dimensions and can be obtained for a wide class of Markov processes. Another 
consequence is the form of the limit - 5' (R'') -valued Wiener process with a simple temporal structure and 
a complicated spatial one. This behaviour contrasts sharply with the case of critical branching systems 
(cf. more detailed description in Introduction). 

AMS subject classification: primary 60F17, 60G20, secondary 60G15 
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1 Introduction 

In this paper, we consider the following branching particle system with immigration. Particles evolve 
independently in R'' according to a time-homogeneous Markov family {Tit,Px)t>o,xm'^- ^he lifetime of a 
particle distributed exponentially with a parameter > 0. When dying the particle splits according to 
a binary branching law, determined by the generating function 

F{s)^q.s'' + {l~q), q<l/2. (1) 

This branching law is subcritical (i.e. number of particles spawning from one is strictly less then 1). 
Each of the new-born particles undertakes movement according to the Markov family 77 independently 
of the others, branches, and so on. New particles immigrate randomly to the system according to a 
homogeneous Poisson random field in R+ x R'' (i.e. time and space) with the intensity measure H\d+i, 
H > (where Xd+i denotes d -I- 1 dimensional Lebesgue measure). Because of immigration the initial 
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distribution of particles does not affect the system in long term. For the sake of simplicity, we choose 
it to be a Poisson random field in R'' with intensity L\d, L > 0. All random objects, the evolution of 
particles, the immigration and the initial distribution are (conditionally) independent. 
The evolution of the system is described by (and in fact can be regarded to be identical with) the empirical 
(measure-valued) process {Nt)t>o, where Nt{A) denotes the number of particles in the set ^ C R'' at 
time t. We define the fiuctuations of the rescaled occupation time process by 

1 

^T{t) ^-rr {Ns- ENs) ds, t > 0, (2) 

Jo 

where T is a scaling parameter which accelerates time (T +oo) and Ft is a proper deterministic 
norming. Xt is a signed-measure-valued process, but it is convenient to regard it as a process in the 
space of tempered distributions S'(R''). The objectives are to find suitable Ft, such that Xt converges 
in law as r ^ +oo to a non-trivial limit and to identify this limit. 

We will discuss some of the related work on the fluctuations of the rescaled occupation time first as it will 
make it easier to understand our result. The series of papers [5, 6, 7, 8, 16, 17, 18, 19] is devoted to the 
study of systems with particles moving according to a symmetric a-stable Levy motion and with critical 
branching (such systems will be referred to as the critical systems, contrary to the subcritical systems of 
this paper). The results therein split roughly into three classes depending on the dimension of the state 
space R"^ 

• "low dimensions" - the system suffers local extinction. The direct study of the fluctuations of the 
rescaled occupation time does not make sense (except from the systems with immigration in [18]) 

• "intermediate dimensions" - the limit has a simple spatial structure (Lebesgue measure) and a 
complicated temporal one (with long range dependence property). 

• "large dimensions" - the limit has a complicated spatial structure (5'(R'*)-valued random field) and 
a simple temporal one (process with independent increments). 

We study the fluctuations of the rescaled occupation time process for systems with subcritical branching, 
which has never been done before. The main result is the functional limit theorem contained in Theorem 
2.1. The functional setting makes the result more interesting and much harder to prove. Moreover, we 
emphasise that the class of systems covered by this theorem is very broad as the restrictions imposed on 
the Markov family ir)t,Vx)t>o^xeR'' mild and natural - cf. Remark 2.1 (this is contrary to the critical 
systems where the proofs rely heavily on the fine properties of a Q-stable Levy motion). In order to 
apprehend the theorem we turn our attention to three aspects of the result. 

Firstly, the theorem is a "classical" functional central limit theorem as the normalising factor is Ft ~ T^^^ 
and the limit is Gaussian, namely a Wiener process. Therefore the temporal structure of the limit is 
simple - the increments of the process are independent. This contrasts sharply with the spatial structure, 
which is a 5'(R'')-valued random field of the form depending on the properties of the Markov family 

Secondly, the subcriticality of the branching law is crucial for the long-term behaviour of the system. The 
limit if of the same nature in all dimensions, making this case much different from the critical systems, 
where the phenomenon of "phase transition" between "intermediate" and "large dimensions" is observed. 
The main reason for this are, roughly speaking, the properties of the movement of particles (recurrence 
vs. transience). In contrast, in the subcritical system the life-span of the family descending from one 
particle is short (its tail decays exponentially) hence the properties of the movement plays much smaller 
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role. Moreover a particle hardly ever visits the same site multiple times which explains the similarity 
of the result to the one for the critical systems in "large dimensions" (i.e. the transient case). It also 
sheds some light on the origin of the temporal and spatial structures. If we consider two disjoint intervals 
which are far away, it is very likely that the increments of the occupation time are contributed by distinct 
families. This results in independent increments. On the other hand the life-span of a family is too short 
to "even out the grains in" which, in turn, gives rise to the complicated spatial structure. We stress that 
the subcriticality of the branching law influences the limit much more then the immigration. The results 
for an analogous immigration system but with critical branching [18] are much different and adhere to 
the scheme for critical systems presented above. See Remark 2.1 for further explanation. 
Finally, notice that the systems considered in the paper do not suffer from the local extinction in "low 
dimensions" (due to immigration each set of positive measure is being populated in arbitrarily large 
times). As it was mentioned already, the limit is of the same nature as in the other dimensions cases. 
It is interesting to compare this results with [11]. "High density" technique applied there enables to 
study the occupation time fluctuations in "low dimensions", which led to a similar limits as in the case of 
"intermediate dimensions". 

To make our paper more comprehensive we present also two illustrative examples in the results section. 
The example one presents perhaps the most important application of Theorem 2.1 to the system of 
particles moving according to a Levy motion. It can be regarded as a subcritical counterpart of the 
critical systems considered in earlier papers (it should be stressed however, that we admit much larger 
class of processes compared to a symmetric a-stable motion previously considered). The resemblance to 
the "large dimensions" case is even more perceptible here - cf. Remark 2.5. In the second example we 
consider a system with particles moving according to the Ornstein-Uhlenbeck process. It is intriguing 
because of the competition of particles attraction towards the origin (caused by the Ornstein-Uhlenbeck 
process) and their disappearance (caused by the subcriticality of the branching law) - cf. Remark 2.7 
Recently, occupation time processes have been intensively. Further to the results mentioned previously 
[9, 10[ presents results for systems with inhomogeneous starting distributions. One should also mention 
[2, 3[ where similar problems are considered in discrete setting (lattice if ). Interesting results were also 
obtained for superprocesses for example [13[ and [12]. In [12] the authors consider a model very similar to 
ours, namely subcritical superprocess with immigration. They only examine the spacial structure (which 
is technically much easier) obtaining a similar result of Gaussian random fleld. One should also mention 
[14] which was pioneering papers in the fleld of systems with immigration. 

The proof technique is similar to the one from previous papers of Bojdecki et al. However, the subcritical 
case required developing new equations and dealing with a general Markov family. This required some 
of the technical arguments to be reflned, as the fluctuations of the occupation time of systems with 
subcritical branching was studied for the flrst time. 

The paper is organized as follows. In Section 2 we present assumptions and the general theorem (i.e. 
Theorem 2.1). Next we give the examples mentioned above. Finally, Theorem 2.1 in Section 3 and 4. 

2 Results 
2.1 Notation 

Before presenting the results announced in Introduction we clear out a few technical points. 5'(R'') is 
a space of tempered distributions i.e. a nuclear space dual to the Schwartz space of rapidly decreasing 
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functions The duality will be denoted by {•, •). 

By (Tt)t>o, A we will denote, respectively, the semigroup and the infinitesimal operator corresponding to 
the Markov family ivt,^x)t>o.xei!i'' presented in Introduction. Sometimes instead of writing ¥.^f{r]t) we 
write E/(j7f ). 

For brevity of notation we also denote the semigroup 

%'^f{x) := e-'^'Ttfix) (3) 
and the potential operator corresponding to it 

/ + 00 
T«/(a;)dt. (4) 



^0 

In the whole paper 

Q = V{l-2q), (5) 

which intuitively denotes "intensity of dying" - recall that V is the intensity of branching and 2q is the 
expected number of particles spawning from one particle. Clearly, subcriticality of the branching law 
implies Q > 0. 

Three kinds of convergence are used. The convergence of finite-dimensional distributions is denoted by 
^fdd- For a continuous, 5'(R'')-valued process X = (Xt)t>o and any r > one can define an 5'(R'''''^)- 
valued random variable 

= ^"{Xt, dt, (6) 

If for any r > X„ — > X in distribution, we say that the convergence in the space-time sense holds 
and denote this fact by — >i. Finally, we consider the functional weak convergence denoted by X„ — >c 
X. It holds if for any r > processes Xn — (Xri{t))teio,T] converge to X = (X(t))tg[o,T] weakly in 
C([0, r], 5'(R'*)) (in the sequel without loss of generality we assume r = 1). It is known that and 
— >/dd do not imply each other, but either of them together with tightness implies -^c- Conversely, —^c 
implies both -^i , ^fdd- 

By c,Ci, . . . ,C,Ci, . . . we will denote generic constants. 

2.2 General case 

Firstly, we present the restrictions imposed on the Markov family {r]t,'^x)t>o,xe'R'i- Not only are they mild 
and quite natural but also are easy to check in a concrete cases (see Section 2.3). First denote quadratic 
forms 

ri(v3) := / (ipix)u'^ip{x)) dx, if, e siR"), (7) 

T2{^) / f/« [T«<^(.)T«f/«<^(.)l {x)dsdx, ^ e SiR"), (8) 

jRd Jo L 

also, slightly abusing notation, we will denote by Ti and T2 the bilinear forms corresponding to them. 
Assumptions 1 

(Al) We assume that the Markov family {rit,Px)t>o,xm'^ almost uniformly stochastically continuous 
i.e. 

Vn sup P4ris,B{x,e)) ^ 1, as s ^ 0, (9) 

X G { — n,n) 

where B{x,e) denotes a ball of radius e with the center in x. Additionally, we assume that for any 
X trajectories of process are almost surely bounded on any finite interval. 
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(Al') Instead of (Al) one can assume stronger but a more natural condition as follows. We assume that 
the Markov family {vt,^x)t>o,xm'^ uniformly stochastically continuous i.e. 

supP^(7?«,B(x,e)) ^ 1, ass^O, (10) 

X 

where -6(2:, e) denotes a ball of radius e with the center in x. 
(A2) Denote by Da the domain of the infinitesimal operator A. We assume 

5(R'*) C Da- (11) 

(A3) Let (fi G 5(R''). We assume that the semigroup (T^)t>o given by 

T,^f{x) = E.exp |-^' <^(77.)ds| fint), (12) 

is a Feller semigroup for any tp. 
(A4) For any ^ £ 5(M'*) 

Ti(^) < +00, T2{ifi) < +00. (13) 



(A5) For any (p € S(R'') 



rp3/2 



[ T^p{x)d.x 0. (14) 



Assumptions 2 

(A6) There exists e > such that for any p 

[ Tt'^p{x)dx<c{lAt-'^-'') . (15) 

(A7) There exists e > such that for any (f and for all h, I 

' T« fT«^(.)7l'^^(-)l Wdx < c (1 A t-'-^) . (16) 
Now we are ready to formulate the theorem which is the main result of the paper. 

Theorem 2.1. Let Xt be the rescaled occupation time fluctuations process given by (2). Assume that 
Ft = T^I'^ and assumptions (A1)-(A5) are fulfilled. Then 

Xr-^iX, and Xr^fddX, (17) 

where X is a generalized S' {R'') -valued Wiener process with covariance functional 

Cov{{Xt,pi) , {Xs,p2)) ^H{sAt) {Tl{(pi,ip2) + VqT2{ipi,ip2)) . 

if, additionally, assumptions (A6)-(A7) are fulfilled then 

Xt ->c X. (18) 

Remark 2.1. Assumptions (A1),(A2),(A3) are typical technical restrictions when dealing with Markov 
processes. We stress that they are mild and fulfilled easily by any "well-behaving" Markov process. 
Condition ( A4) is natural, as it states only that the limits is well defined (if it is not fulfilled normalization 
larger then Ft = T^^^ is required). Finally to analyse (A5) let us notice that 

/ %'^lA{x)dxds, (19) 

JR'' 



5 



is the average number of particles in set A for the system starting from the null measure. Intuitively, 
the aim of (A5) is to "prevent gathering infinite number of particles" in any set. Taking this explanation, 
(A5) seems to be too strong as T~^~'^ for any e > should be sufficient. 

To sum up we state questions which raise naturally for further investigation. Firstly, the gap between 
(A5) and T^^^'- is somehow unpleasant. A natural conjecture is that Theorem 2.1 holds also with this 
weaker condition. Another, less likely in the author's opinion, possibility is that the gap can be explained 
in probabilistic manner. Any result in this field would possibly give Theorem 2.1 even more elegant form. 
Secondly, assumptions (A3),(A4),(A5) impose a certain regime of behaviour on the system, in which the 
subcriticality suppress the contribution of the motion to the limit. By relaxing them the contribution of 
the motion "increases". Rough calculations suggest that this in turn results in an increase of the norming 
factor Ft and in the limit with a complicated temporal part. However, with the motion playing larger 
rdle, this case it is not likely to be captured as generally and elegantly as in Theorem 2.1. The reader is 
also referred to Remark 2.7 for more detailed explanation in a particular example. 

Remark 2.2. Assumptions (A4),(A5) are clearly technical. It is not obvious whether they are necessary. 
This question have not received enough attention yet, as the main goal of this paper was to identify the 
limits. Finding necessary conditions for tightness seems not to be easy, however. 

Remark 2.3. As it was mentioned in Introduction the limit is a 5'(R'^)-valued Wiener process with 
a simple time structure and a complicated temporal one in all dimensions. This result resembles the 
result for the system with critical branching in large dimensions. The main reason of this is a short 
(exponentially-tailed) life-span of a family descending from one particle. On the one hand it leads to 
independent increments in the limit (as there are no "related" particles in long term). On the other hand 
the movement is "not strong enough" to smooth out the space structure. 

Another remarkable, yet not such unexpected, feature is that the limit can be obtained for "low di- 
mensions". Due to immigration the system no longer suffers from local extinction and the limit can be 
obtained without special techniques, like high density limits of [11]. 

Remark 2.4. Analogous systems but with critical branching were studied in [18]. The results there 
are much different from Theorem 2.1 and adhere to the scheme observed for other critical systems (as 
described in Introduction). This proves clearly that the subcriticality influences systems much more then 
the immigration, at least with respect to the limit behaviour of the fluctuations of the occupation time. 

2.3 Examples 

The theorem in the previous section is quite abstract. Now we will present two illustrative examples. 

Levy motion Recall the description of the system A'^ from Introduction. In this example the movement 
of particles will be given by a Levy process - we keep the notation, by {r]t)t we denote the Levy motion 
starting from 0. Its characteristic function is 




(20) 



where ^' is the Levy-Khinchine exponent i.e. 




(21) 
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where a G (drift term), K is non-negative defined n x n matrix (covariance of the Gaussian part) and 
/u is a (spectral) measure obeying condition J^-^^gyi^^ A l)da; < +00. 

Let us now check that the Levy motion fulfils Assumptions 1. It is a space homogeneous process hence to 
check (Al') it suffices to show that rjt -^^ which follows directly from the characteristic function. (A2) 
is slightly more difficult - let us take (f € 5(R''); one can check that 



t 

which implies 

i'i'{z)(p{z). 



Ttlfi{z) - (p{z) 



t 

Hence </3 € Da- 

We skip the proof of (A3) which is as an easy consequence of the space homogeneity and Lebesgue's 
dominated convergence theorem. 

Recall that A is an invariant measure of the Levy motion. We have 



(22) 



It can be checked that U'^Lp[z) — (p{z){Q — ^{z)) ^ . Applying the Fourier transform to (7) we obtain 
T2 can be treated in a similar way 

The real part of ^ is non-positive hence clearly both Ti{(p) and T2{ip) are finite, therefore (A4) holds. 
The assumption (A5) follows easily from calculations below 



-QT 



Finally (A6)-(A7), can be proved in the same way. Utilising Theorem 2.1 we obtain 

Theorem 2.2. Let Xt he the occupation time fluctuation process given by (2) for a system of particles 
moving according to a Levy motion. Then 

Xt -^c X, asT ^ +00, 

where X is a generalized S' {R''') -valued Wiener process with covariance functional 

Cov{{X,,,fi,),{Xt,^2)) = 

(^"^)f(^I(Q^+(Q^w)^^^^^^'^ (25) 
Remark 2.5. Formally the result resembles the result for the critical branching systems in "large dimen- 
sions". Indeed, by converging with branching law to a critical one (i.e. q ^ 0) and decreasing intensity 
of immigration (i.e. H ^ 0) appropriately in the rhs of the expression above, one gets 

^ ^) (2^)^ L + (2^) 

which exactly the limit in theorem [6, Theorem 2.1] (with ^{z) — —z" for the symmetric a-stable Levy 
motion considered there). The question if this convergence has any probabilistic meaning is natural but 
has not been addressed yet. 
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Remark 2.6. Denote the spacial part of the limit by Y. It easy to notice that it is a homogeneous 
(generalized) Gaussian random field. The measure 

M(d^) := (ttA^ + TTT-^XTTTal (27) 



'-*(2) {Q-^{z)f 

is called the spectral measure of Y. It is well-known (see e.g. [15, Preposition 1]) that Y is "classical" 
i.e. is function-valued random field if and only if its spectral measure is finite. For the system considered 
in this section this translates to the condition 

/ TT-^-rr^dz < +<x. (28) 

In the most important case when the particles move according to the symmetric a-stable Levy motion 
the condition is true if and only if d = 1 and a G (1,2]. 

Ornstein-Uhlenbeck process In this example the movement of particles is governed by the Ornstein-| 
Uhlenbeck process. The Ornstein-Uhlenbeck process is the solution of stochastic equation 

dXt ^ -eXtdt + adWt, 6) > 0,(7/0. (29) 

Its semigroup is given by 

Ttfix) = (5„„(t)/) (xe-"') (30) 

where 5 is the semigroup of Wiener process and ou{t) = (1 - e~^**)/(26l), OU{t) = (e^** - l)/(26'). 
Assumptions (A1),(A2) i (A3) can be checked easily from the following representation of the Ornstein- 
Uhlenbeck process 

Xf — xe H ■=W (e —lie 

V2e 



Recall definition 5. We assume also that Q > 9 (this assumption is crucial and will be explained later in 
Remark 2.7). We have 

/ Tt'^f{x)dx = e-es-"'* / f{x)dx. 
Using this equation (A4)-(A6) can be easily verified. 

Using the Fourier transform we can calculate Ti and T2 in more explicit form. Namely, by (30) and the 
fact that the Lebesgue measure is an invariant measure of S we have 

Ti{ip) = I U'^ (ip{x)U'^ <p{x)) dx = ^ 1 <fi{x)U'^ip{x)dx. 
Using the Fourier transform we get 

Similar calculations for T2 give 

+^ /"+°° ^_(Q_9){2.+.) f ^(e«^^)e-20^(»)l^lV(e«*z)e-°^("'l^''°^lV(e«(^+")z)dzdnds. 
[2-k) Q Jq Jo Jm<' 

(Recall that quadratic forms Ti and T2 induce corresponding bilinear forms). Assumptions (A6), (A7) 
can be easily verified this entitles us to use Theorem 2.1. 
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Theorem 2.3. Let Xt be the rescaled occupation time fluctuation process given by (2) for a system of 
particles moving according to the Ornstein-Uhlenbeck process. Assume that Q > 9. Then 

Xt X, asT ^ +oo, 

where X is a generalized S' (R'') -valued Wiener process with covariance functional 

C0V{{X,,^^) ,{Xt,ip2)) = («^*)^(^ (ri(v3l,</92) + l/<?T2(<y3l,<^2)) , </9l , ^2 G 5 (r'') . 

Remark 2.7. This example is interesting because we can observe struggle of two antagonistic forces. One 
is the "exponential attraction" of particles from the whole space to the proximity of by the Ornstein- 
Uhlenbeck process the other is dying out of particles because of the subcriticalty of the branching law. 
More precisely, denote = iB(o.r) then 

/ T«<^(^)dx = e-«* / (5„„(i)</p) [xe-'')Ax = e-(«-'')*|B(0, r)\. 

is the average number of particles in the ball -6(0, r) for the subcritical system (without immigration) 
starting out from the homogeneous Poisson field. The condition Q > 6 can now be easily interpreted - 
the subcriticality is "strong enough" to prevent gathering of particles (near 0). 

This observation raises a natural question what happen when Q = 9 i.e. when the forces are in the perfect 
balance. Rough calculations suggest that norming factor is greater ( Ft = T) and the properties of the 
motion affect the temporal part of the limit (it is not longer process with independent increments). 

3 Proofs 
3.1 Scheme 

To make the proof clearer we present a general scheme here and defer details to separated sections. 
Although the processes Xt are signed-measure-valued it is convenient to regard them as processes with 
values in 5'(K'*). In this space one may employ a space-time method introduced by [4] which together 
with Mitoma's theorem constitute a powerful technique in proving weak, functional convergence. 

Convergence Prom now on we will denote by Xt (r = 1) a space-time variable (recall (6) with r = 1) 
defined for Xt . To prove convergence of Xt we will use the Laplace functional 

Lt(4') =Ee-<-*^'*), $ G 5(R''+'),<E> > 0. (31) 

Por the limit process X denote 

L($) =Ee-<-^'*), $ G 5(R''+'),$ > 0. 

Once we establish convergence 

Lt{'^) L{^), asT^ +00, V4.g5(Rd+i),,j>o. (32) 

we will obtain week convergence Xt X and consequently Xt X. Two technical remark should be 
made here. We consider only non-negative $. The procedure how to extend the convergence to any $ is 
explained in [5, Section 3.2]. Another issue is the fact that {Xt, "J"^ is not non-negative (which is a usual 
condition to use the Laplace transform). The usage of the Laplace transform in this paper is justified by 
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the special (Gaussian) form of the limit. For more detailed explanation one can check also [5, Section 
3.2]. 

As explained in [7] due to the special form of the Laplace transform convergence (32) implies also finite- 
dimensional convergence. 

Detailed calculations for this part of scheme will be conducted in Section 3.3 and Section 4.1. 

Tightness Using additional assumptions the tightness can be proved utilizing the Mitoma theorem 
[20], which states that tightness of {Xt}j- with trajectories in C([0, 1],5'(R'^)) is equivalent to tightness 
of {Xt,4>), in C([0,r],M) for every 4> G 5(R''). We adopt a technique introduced in [6]. Recall a classical 
criterion [1, Theorem 12.3], i.e. a process (j) e 5(R'') is tight if for all t,s>0 

E({XTit), ^) , {Xt{s), ^))* < C{t - sf. (33) 

Following the scheme in [6] we define a sequence {ij)n)n in 5(R), and x»i(^) = /J V'n('S)ds in a such way 
that 

0<Xn <![»,*]■ (34) 

Denote = ® i/i^ . We have 

lim {XT,$n) = {Xrit),^) - {Xt{s),^) 



n — ' + 00 



thus by Fatou's lemma and the definition of V'n we will obtain (33) if we prove (C is a constant independent 
of n and T) that 

E^Xt,*™)* < C{t-sf. 

From now on we fix an arbitrary n and denote "I> :— $,1 and x X"- By properties of the Laplace 
transform we have 



Eexp ( ~e{XT,^ 

9=0 



Hence the proof of tightness will be completed if we show 



Ee^p i-e ( XT,'i>) ] <Cit-sf. (35) 



9=0 

Further calculations are deferred to Sections 3.4 and 4.2. 

3.2 Auxiliary facts and one-particle equation 

In this Section we will prove an equation for one particle which will play a key role in the rest of the 
proof. Before that we recall a general Feynman-Kac formula which is crucial for its proof. 
Let A be a (unbounded) linear operator with domain Da- We define a problem 

^w{t,x) ^ Aw{t,x) + c{x)w{t,x), ^^^^ 
w(0,x) = f{x), 

where w{-,t),f £ Da- 

Proposition 3.1 (Feynman-Kac formula). Let (Xt,Pi) be a uniformly stochastically continuous Markov 
family (cf- assumption (Al)) with values in with infinitesimal operator A- Assume also then c is a 
uniformly continuous and hounded. Then 



■u: 



{t,x)=E^expl c{Xs)ds} f{Xt), t>0,xeE 
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is a solution of (36). It is the only solution in the class of functions w such that sup^ \x{x,t)\ < e"', Vt 
for Q e M. 

Recall that F denotes generating function of the branching law (1). We define G(s) = _F(1 — s) — (1 — s) 
so in our case 

G(s) = qs"^ + (1 - 2q)s. (37) 
Behavior of the system starting off from a single particle at x is described by the function 

(a;,r,t) = 1 -Eexpj-^ (iV,"^, * (•, r + s)) dsj , * > 0, (38) 

where denotes the empirical measure of the particle system with the initial condition A^q = Sx- More 
precisely A^^ is a system in which particles evolve according to the dynamics described in Introduction 
but without immigration. 
The lemma gives the announced equation 

Lemma 3.1. Assume that ^ > and assumptions (A1)-(A3) are fulfilled then 

«*e[0,l], (39) 

and Vit satisfies equation 

{x, r,t)^ I Tt-s [* (•, r + t - s) (1 - (•, r + f - s, s)) - VG (u* (•, r- + t - s, s))] (x) ds. (40) 
Jo 

Formally, this equation is the same as [5, (3.22)]. We have to refine the proof as in this paper we 
consider a more general case hence 

Proof. (39) follows directly from the definition (38). 
Now we proceed to the proof of (40). Denote 

Wijiix, r, t) := 1 — vmix, r, t) 

In the first step we expand (A3) to a slightly more general semigroup. For ^ e 5(R''^^), r > define 

T*'7(x) = Eexp |- ^' (ry:, * (-, r + s)) dsj /(r,f ) 
We claim that T* ' is also Feller. 

Define ^n{x,t) = X^fc^i ^'''(^^j *fe)l[tfc_i,ti,) (i) where ~ tk/n. Inductive argument (with respect to n) 
implies easily that T'^""'~ f{x) is continuous when / is continuous. Indeed one can write 

T*-7(a:) = exp |- *(7)s,r- + ti)| exp |- *„(r?,,r + s)ds| . (41) 

Using the Markov property we have 

T*-7(a;) =E.exp|- r*(,,.,r + to|T*-'-+'V(^tJ =T*'-''^+'-)T*?;'-+*V(^). (42) 



By induction we can assume that T^_^'_^^^'^ f{x) is continuous and by the Feller property of 7^*( '*i+'') 
(assumption (A3)) we get asserted claim. It is obvious that 7^*"''^/ —> T*'*"/ uniformly hence T^'^f is 
continuous. 

In the next step we will prove that this fact implies the continuity of wxa. Denote by T the space of 
ancestor trees i.e. a space of binary trees with nodes and leafs labeled by the splitting and death times 
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of particles respectively. The splitting dynamics described in Introduction induces a probability measure 
on them - v. Let us notice that in our case the trees are finite almost surely. For a given ancestor tree r 
we can calculate Wt,* given by 



:= Eexp |- ^ (M^"', *(•, r + s)) dsj , 



where is the branching particle system with the branching dynamics encoded by r. Let |r| denotes 
the height of r. By induction with respect to the height of the tree we can prove that Mt,* is continuous. 
For trees of height 1 it is obvious from the continuity of z. Let r be a tree such that n = |r|. Removing 
the root splits r into two sub-trees n, r^. By t\ we denote the label of the root i.e. the time of the first 
branching (or death if the root is a leaf). If ti > i the continuity if obvious hence we are remain only 
with the situation when t\ < t. One can write 

■WT,9{x,r,t) = Eexp |— J {rj^, + s)) ds| (wri,*(%i,r + U,t — U)wT2,^{rhi,r + ti,t — ti)) . 

Now continuity of Wt,<!/{-, r, t) follows from the induction hypothesis and Feller property of T*'''. Further, 
it can be easily proved that Wt,* is in fact continuous as a function of three variables. This property infer 
continuity of lo* which is justified by the formula 

■w<s,{x,r,t) = j Wr,<s{x,r,t)v{AT). 

and Lebesgue's dominated convergence theorem (recall < Wr,* < 1 and is a probability measure). 
Conditioning on the time of the first branching we get 



w 



+ V j e^'^'Eexp^-^ '<il{r]l,r + u)Av^ F {w{ru,r + s,t- s)) . 



Let us introduce functions h i k 



h{x,r,t)~e ^*Eexp|-^ * (r?^' , r + s) ds 



- ^ ^{riu,r + u)du^ F [w {r,t,r + t, s)) . 



ks {x,r,t) ~ e ^'Eexp 
Now w can be written as 

w {x,r,t) = h {x,r,t) + V ks {x,r,t — s) ds. (43) 
Jo 

The crucial step of the proof is application of Feynman-Kac formula. Assume for a moment that the 
Markov family fulfils (Al') instead of (Al). Let 9 £ 5(R''+^) and define 



l0{x, r, f) = E exp 



Assumptions (AT), (A2) assert that we can use Proposition 3.1 (one have to prove that O belongs to the 
domain of the infinitesimal operator of Markov family t — > {Xf,r + t) - we skip this simple step) hence 
le{x,r,t) is solution of 



§-Je{x,r,t) = (A^ + ^ - 1'(a;,r-)) le{x,r,t), 
l0{x,r,O) = 6(2;,?-). 



12 



Let us denote 

keix,r,t) = e~^''le{x,r,t). (44) 

Direct computations yield 

j-^ke{x,r,t) = (A„ + ^) ke{x,r,t) - {■^{x,r) + V)ke{x,r,t), 
ke{x,r,0) — Q{x,r). 

This is an evolution equation which has an integral form 



ke{x,r,t) =Tte{x,r + t) - f Tt^s[{'l'{-..r + t-s) + V)ke{-,r + t-s,s)]{x)ds. 

Jo 



(45) 



Now define r„ = inf{t : \rjt \ > n} and processes r}"'^ ~ Vr^At- Clearly they are Markov and each of them 
fulfils (Al) and (A2). We know so far that 

k^{x,r,t)^Tt"e{x,r + t)~ f Tr_sl{'^{-,r + t^s) + V)k'^{-,r + t-s,s)]{x)ds, (46) 

Jo 

where T",kQ denote respectively semigroup and equation (44) defined for Markov process rj^'^. It easy 
to show that k"^ ^ k (point-wise) and consequently using Lebesgue dominated convergence theorem show 
that (45) is fulfilled for any Markov family satisfying (Al). 

Clearly F{w{-, •, s)) is continuous hence there exists a sequence {On)n, Qn G 5(R'*) convergent uniformly 
to F{w{-,-,s)). Applying this to definition (44) we obtain point-wise convergence 

fce„ {x, r, t) — > ks{x, r, t). 

Now we use dominated Lebesgue's convergence theorem {ke„ < sup 6„ < c) to the right side of (45) 

ksix,r,t) ^TtF{w{x,r + t,s))- [ Zs [('i'{-,r + t - s) + V) ks(-,r + t - s, s)] {x)ds. (47) 



Analogously 

h{x,r,t) = l- [ Tt-s[{<i'{-,r + t- s) + V)h{-,r + t- s,s)]{x)ds. (48) 
Jo 

We put the obtained equations in (43) 

w{x,r,t) = l- [ Tt-s[{'l'{-,r + t-s) + V)h{-,r + t-s,s)]{x)ds+ 
Jo 

V [ Tt^sF(w{x,r + t-s,s))ds- 
Jo 

rt rt — s 

V I I Tt-s-u [(*(•, r + t- s-u) + V)ks{-,r + t- s-u,u)] {x)duds. 
Jo Jo 

We substitute u ^ u — s and change the order of integration 



w{x,r,t) = l- / Tt-s[('^{-,r + t~s) + V)h{-,r + t-s,s)]{x)ds+ 

V [ Tt-sF{w{x,r + t- s,s))ds- 



V / ks{-,r + t — u,u — s)ds 



{x)du. 



Tt_„(*(-,r + t-u) + V) 
Finally we apply (43) to the second and fourth term 

w{x,r,t) = 1 - I Tt-s[{-^{-,r + t-s) + V)w{-,r + t-s,s)]{x)ds+ 
Jo 

V [ Tt-sF{w{x,r + t- s,s))ds. 
Jo 



13 



Recall that 1 — w = Finally trivial computations yield asserted (40). 

□ 

We consider the case of subcritical branching {q < 1/2) in (1). Recall equation (37) and Q — V{l~2q), 
putting this to equation (40) gives 

v^{x,r,t)= [ Tt-s['i'{-,r + t- s){l-v^{-,r + t- s,s))- 
Jo 

Qd* (•, r + t — s, s) — Vqvxi (•, r + t — s, s)^] (a::) ds. (49) 
vi is quite cumbersome to deal with hence we will approximate it with -C* defined in the following way 

V9{x,r,t) ^ [ Tt^^<l'(-,r + t- s)ds, * G 5(R''+^), x £ R'', r, t > 0. (50) 
Jo 

It can be easily checked that this function fulfills the equation 

V9{x,r,t)=f Tt-s['i'{-,r + t-s)-Qi}^{-,r + t-s,s)]{x)ds. (51) 







Intuitively -5* was obtained by dropping quadratic terms in (49) which "do not play role" when 'I' is small. 
The quality of the approximation is expressed in terms of function u 

:= — (52) 

We have 

Lemma 3.2. Let ^ > 0, then satisfies the equation 

u^{x,r,t)= [ Tt1^['i'{-,r + t-s)v^{-,r + t-s,s) + Vqvl{-,r + t-s,s)]ds. (53) 
Jo 

Proof. Subtracting equations (49) and (51) we obtain 

UiS!{x,r,t) = / Tt-s (5ii*(-, r + t ~ s,s) + 'i'{-,r + t ~ s)u*(-, r + t — s,s) + Vqv^{-,r + t~s,s)] ds. 
Jo 

(54) 

Although we do not know solution of (49) we may treat as a known function. It is easy to check 
that (53) solves (54). Standard application of the Banach contraction principle proves that that it is 
unique. □ 

Notation For now on we fix non-negative $ and prove convergences announced in the scheme in 
Section 3.1. To make the proof shorter we will consider $ of a special form 

^{x, s) = ip{x)tp(s), ip e 5(R'*), tp e S{R),ip>0,tp >0. (55) 



We also denote 



We write 



VT (x) = -^(fi (x) , x(s) = ^ V(^)d«, XT = x(j^] ■ (56) 



^{x,s) = (p{x)x{s), 
^t{x,s) ^ —'i' (x,-^) =ipT{x)xTis). (57) 



Ft \ T 

note that ^ and 5't are positive functions. In the sequel, we also write 



VT{x,r,t) = Vi!j^{x,r,t) and vt{x) = vt{x,0,T), (58) 
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and 

It is obvious now that ?i* > which together with equations (52) and (50) implies 

Ft 

We will also use the following simple estimation 



(J 

UT{x,r,t) < — 2 



(59) 



(60) 



Fix ^ and denote 

v{e) = ve9 (61) 

In the sequel we will need derivatives of v with respect to 6. It is easy to calculate by (49) that (we omit 
arguments and integration variables) 



When 9 = then 



v'{6)= f Tt-s [<f{i~v(e))~eW{e)-Qv'{e)~2Vqv{9)v'{e)] 

Jo 

v'{0) = pTt-s [*-Qi/(0)] . 
Jo 



(62) 



(63) 



It is easy to notice that it is the same equation as (51) hence v = v'{0) (note that the above calculation 
is not quite rigorous as one have to justify differentiation under integral in (62)). 



3.3 Laplace transform 

In this section we calculate the Laplace transform of the space-time variable Xt. Let us recall that the 
initial distribution is given by a Poisson random field with intensity LA, L > and the immigration is 
determined by a Poisson random field on R+ x R'* with intensity H {X^ X) , H > 0. We can split the 
system A*" into two independent parts 

A^t =iV° + iV/""", 

where A''" consists of particles being in the system at time t = and their offspring while n'"^"^ is the 
immigration part with particles which appeared in the system after t = and their descendants. 
The first step is calculating the Laplace transform of the space-time variable corresponding to rescaled 
occupation time process Y. It is easy to check that 

{NTs,<f{-,s))ds = / {Ns,<fT{-,s))ds 



Denote (recall (57) for the relation between ^ and <&) 

Kt{^) = Eexp (Yt, = Eexp ( - 

We can write 

AV($) = Eexp |- f {N°, *t(-, s)) ds J-Eexp 



7V/'"'",*T(-,s))ds 



(64) 



(65) 



(66) 
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Firstly evaluate the term with N'''^^ . Conditioning with respect to Imm, using independence of evolution 
of particles {branching Markov property) and equation (38) we obtain 

E(^exp|-^^(^Arf™,*T(-,s))ds| Imm^ = 

n Ee^pi^-J^ (iV,^_t,*T(-,s))ds| = n {l-v^Ti^,i^T~t)), (67) 

{t,x)fzZmm it,x)^Imm 

where Imm is a (random) set such that X^^^ x)<Umm. ^(t,^) ~ ^n^n^os- and 5(t,a;) corresponds to a particles 
which immigrate to the system at time t to location x. Hence we have 

Eexpj-^ <^Ar/'"",*T(-,s)^ds| =Eexp {(/mm, log(l-«T(-,*,r -★))}, 

where denote integration with respect to space and time, respectively. Taking into account distribution 
of Imm we obtain 

rT , ,1 r rT 



Eexpj-^ (:^iV/'"'",1'T(-,s)^ds| = exp|-_ff^ j (a;, T - t, t)da:dt j 
ad can be treated similarly - we have 

{N°,^T{-,s))ds^ =expi^-L j u*^(a:,0,r)dx| . (68) 



The first term is easier and can be treated similarly - we have 
Eexp 



Finally we have 



Kt{^) =expl-H f f v^,^{x,T -t,t)dxdt~ L f v^,^{x,Q,T)dx 



(69) 



By the properties of the Laplace transform we have (recall also that i't(O) — vt - see (63) and v{Q) — 0) 



E(y.,<.)=A 



KT{9^) = ~Hf f v<i^{x,T -t,t)dxdt - L I iq,^{x,0,T)dx. 

?=0 Jo JR'^ 



Now we can calculate the Laplace transform of Xt- Using definition of ut (52), simple fact that Xt = 
Yt — EFt we obtain 

rT 



(70) 



Lt($) = Eexp j-fJ^Xr,*^ I =cxp|iy UT(a;, 0, r)ds:| exp ji? ^ J ut{x,T - t,t)dxdt 
Now the task is to show limit of (70). Using (53) one obtains 

Eexp (Xt, '!>)}= exp{L (Ai(r) + A2{T)) + H {A-i{T) + Ai{T))} , (71) 

where ^ 

Ai{T)= [ [ T^_^['i'T{-,T ^ s)vt(:,T ~ s,s)]{x)dsdx, (72) 
Jm.'^ Jo 

A2{T)^Vqf f T^_^Vt{x,T - s,s)dsdx, (73) 
Jr'^ Jo 

MT)= f f f T,'^^A'^T{-,T-s)vT{-,T~s,s)]{x)dsdtdx, (74) 
Jw^ Jo Jo 

A4{T) = Vqf f f Tt'^^v'^ix,T - s,s)dsdtdx. (75) 
Jr^ Jo Jo 

The first part of Theorem 2.1 will be proved once we have established 

Ai(r) ^ 0, A2(T) -> 0, asr^+c50. (76) 
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MT)^ [ x(l-s)'/' U'^\ip{-)U'^<p{-)\{x)dxds, asT^+oo. (77) 

Jo JRii L J 

A4{T)^2Vqf x{i~vif r [ U'^ \T,'^^{-)T,'^U'^f{-)]{x)dxdsdvi, asT^+<x. (78) 
Jo Jo Jrii -' 

In the next section we will prove (77), (78). The poofs of (76) are simpler and are left to the reader. 

3.4 Tightness 

Recall that we continue the proof according to the scheme in Section 3.1. First, we will compute the left 
hand side of (35). We adopt the following notation - denote $e,T = d^T and '^b.t ~ = 9ipT ® xt 
related to $e,T by equation (57). Additional parameter will indicate that a particular quantity is 
calculate for $e,T or '^g^r- Hence using (71) we can write (we additionally assume that L = 0, H = 1, 
the proof without this assumptions goes exactly the same lines but is longer) 

EexpS^-(^XT,e^)j =exp{A3{e,T) + A4{d,T)}. (79) 

For sake of consistency we denote 

v{6) ■— v{6){x, r, t) = y {x, r, t), 
Differentiating (49) and evaluating at S = yields (we skip arguments and integration variables) 

v{Q) = 0, 

ru 

v'{Q)= / T„_, [*T"O«'(0)] , 
Jo 

ru 

v"{Q) = / T„_, [-2*T^'(0) - Qv"{e) - 2Vqv'{0f] , 
Jo 

f'U 

w"'(0)= / T„_. [-3*Ti'"(0)-Q7;"'(0)-5Vgi;"(0)i;'(0)] . 
Jo 

These equations can be solved (we skip detailed calculations) 

w'(0)(s:,r,t) = I T^_^^T{x,r + t-s)ds. (80) 
Jo 

v"{0){x,r,t) ^ ~2 [ T,?, [<i'T{x,r + s)v'{0){x,r + t- s,s) + 
Jo 

Vqv'{0){x,r + t- s,sf]ds, (81) 

v'" iO)ix, r,t) = - / T^l^ [3'^T{x,r + t-s)v"{0){x,r + t-s,s) + 
Jo 

5Vqv" {0)(x, r + t — s, s)v'{0)(x, r + t — s, s)] ds. 

Differentiating equations (74) and (75) and evaluating at 6 = one gets (in the last expression we skip 
arguments) 

^3(o,r) = o, a;,(o,t) = o, 

A^\Q,T)^i [ f f T,'^J'(i'T{x,T-s)v^'-'-\0){-,T-s,s)\{x)dsdtdx, i>2. 
Jr'' Jo Jo '- ^ 

A4(o,r) = o, Al(o,T) = o, 
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A2{0,T)^2Vqf [ [ %'^XOy{x,T - s,sfdsdtdx. (82) 
Jm^ Jo Jo 

^' (0, T) = Vq [ r f («(0)'%(0)' + {v{Q)"f) AsAtAx. (83) 

Jp.rf Jo JQ 



exp {Az{e, T) + Ai{e, T)} = 4^(0, T) + ^^(0, T) + 3(A:,'(0, T) + <(0, T))= 



Now we are ready to differentiate (79) 
Now in order to show (35) it suffices to prove 

^r(o,T) <c(^-s)^ Ar(o,T) <c(^-s)^ (84) 

A'i'{Q,T)<c{t~s), A'i'{Q,T)<c{t~s). (85) 
Example computations will be shown in Section 4.2 

4 Calculations 

4.1 Calculations - convergence 

Convergence of A^, Firstly, we replace v with v. Secondly, we calculate the limit for such expression. 
In the end we will prove that the change do not affect the limit 

MT)= f f f T,'lA-^T{-,T-s)vT{-,T-s,s)]ix)AsAtAx. 

jRd Jo Jo 

Using equation (50) and Fubini's theorem we get 

f-T ft 



III r T,'i\<i>T{-,T-s)T^_^-iiT{-,T-u)](x)AuAsAtAx. 
JvA Jo Jo Jo ^ 



10 JO JO 

Using (57) and Fubini's theorem once more one can write 

rT rt rs 



MT)= / / xt{T-s)xt{T-u) %'iAvT{-)T^^^VT{-)\[x)AxAuAsAt 

J rt J n J n Jwd L J 



10 Jo Jo 

Changing variables t Tt, s Ts, u Ts and using (56) we have 

i3(r) -^m^ X(l - ^)X(1 - u) j^^ T«,_^) [p{-)T^(s-u)Vi-)] {x)AxAuAsAt. 
Changing the order of integration and changing u s ~ h one obtains 

A;{T) = x(l ~ s)x(l -' + ^^ [ j^, '^nt-s) [fi-)'^T,M-)] {^)dxAtAhAs. 

Finally changing t ^ t + s we obtain 

MT) = S / V - -^ + h) r W T« Li-)T^,M-)\ ix)dxdtAhAs. 

l^T Jo Jo Jo JvLd L J 

Now it is obvious that 

A^{T) = A3i{T) + A;2{T), (86) 
where , ^ ^ 

M{T) r X(l -^f r I V '^Tt U-)T^M-)] ix)AxAtAhAs. 

Jo Jo Jo JRd L -' 
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Az2{T) f x(l - s) [x(l -s + h)-x{l- s)] ' T« [^{■)T^M-)] {x)dxdtdhds. 



Ft jo 



Recall Ft — T^^" and change integration variables t — + t/T and h h/T 

r-l ,-Ts ,-T{l-s) 



Mi{T)^ ( xi^-sfT I ' / T« L(.)T«^(.)1 (a;)da;dtdMs. 

Jo Jo Jo Jm<' J 

Fubini's and Lebesgue's monotone convergence theorems imply immediately 

Mi{T)^ [ xO--sf [ ^U'^[^{-)U'^^{-)\{x)dxds, asT^+oo. (87) 



Notice also that by assumption (A5) the integral is finite. For a (5 > one can choose e > such that 
suPhs(o,£) lx(l - s + ^) - x(l - s)| < <5 we have 

\M2{T)\ < 5^ r r r W j-Q L(.)r«^^(.)] {x)dxdtdhds+ 

^T JO Jo Jo JR'' L J 

^ r r r ^ f L(.)t«,^(.)1 (xodx-d^ws 

^T JO Jt Jo JR'' L J 



By virtue of convergence (87) we know that the first integral is finite the second can be written as 

^T Jo Jo Jo JW^ 

Changing integration variables h h/T and t t/T and using Ft = T^'^^ we have 

"1 /' + OC /' + 00 



l(,,,)(/i) f l(o,i-s)(i) / ■T't* ['/'(•)^T'i<^(-)] (a;)dxdM/ids = 



1 

^0 Jo Jr'' 



(T„,TB){h)l(0,T{l~s)) {t)T,'i [^(•)T«<^(-)] {x)dxdtdhds 



The integrand converges point-wise to and is dominated by T^'^ |^<^(-)T^''(^(-)j (x) (which by virtue of 
previous argument is integrable) hence by Lebesgue's theorem the integral converges to 0. 5 can be chosen 
arbitrarily small, consequently 

i3i(r)^0. (88) 
The last step is to estimate the difference As{T) — A3{T). Be definition (52) we have 

A3{T)-A3{T)= [ [ [ %'^^['fT{-,T - s)ut{-,T - s,s)]{x)dsdtdx. 
jrA jo Jo 

We can utilize inequality (60) hence 

A-s{T)-As{T) <^ f r r TSs'fT{x,T - s){x)dsdtdx. 

^T JR.'' Jo Jo 

Using notation (57) we have 

AsiT) - AsiT) < H / r f Z'lM^)dsdtdx. 

^T JR.'' Jo Jo 

By assumption (A5) it is straightforward to check that this converges to as T +oo. 
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Convergence of A4 Similarly as for A3 we replace v with v and calculate the limit for such changed 
expression. In the end of the section we will prove that the change do not affect the limit 

A4iT) = Vq f r /V,! [vl{.,T-s,s)] {x)dsdtdx. (89) 

jRd Jo Jo 

Firstly we use equation (50) 

A4T) = Vqf r f r /V.^, [T5,^*T(-,T~i;i)T«,^*T(-,T-t;2)l (a;)d«2d«idsdMx. 
Jmd Jo Jo Jo Jo L J 

Using (57) and Fubini's theorem yield 

l-T pt 



A4T)^Vq / / / Xt{T-vi)xt{T~V2) T,'l\T^'i^^VT{-m'l,,M-)]{^)dxdv2dv 



idsdf. 



10 Jo Jo Jo Jp.^i 
We substitute t Tt, s —* Ts, vi — > Tvi, V2 Tv2 and use (56) 

rp4 /■! /"t /"S /"S 

^ , , , , 

T Jo Jo Jo Jo 



A4T)^Vqj^ /, /, / /, x(l-^i)x(l-^'2) y^^T«,_,) [t«_^)<^(.)T«_^j^(.)] (x)dxdi;2d«idsdt 
Next we change the order of integration and use symmetry 

A4T) = 2Vq^ r r x{l - vi)x{l ^ V2) 
r-r Jo Jo 

f f I [t« v'(-)7't1.-„,)^(-)1 (a:)d^dtdsdi;2dt,i. 

Jvi J s JvA -' 

Let us now substitute V2 ^ vi — s ^ s + vi^ t ^ t + s 



Ai{T) = 2Vq'^ f r x(l - ^'i)x(l - ^'1 + h) 

Jo Jo 

' "'^^T^^ [t^«,^(.)T^'^^+;,)^(-)] (^)dxdtdsdMw. 

Now it is obvious that 

Ai (T) = ^41 (T) + Ai2 (T) , (90) 

where 

? , , , , 

T Jo Jo Jo Jo 



i4i(r) = 2Vqj^ / x(l - ^i)' / / / T^t [t«<^(.)T«^^,,^(.)] {x)dxdtdsdhdv 



T 



A4i(T) = 2Vq^ f r x(l - t'i)(x(l - ^'1 + /i) - X(l - ^'1 + /i)) 
I't Jo Jo 



l — vi pl~ s — vi 



Jo 



T ' 



[T^M-)T?(s+hM-)] (a;)da;dMsdMi;i. 
Recall that Ft = T^/^ and substitute h ^ h/T, t t/T, s ^ s/T 

Aii{T) = 2Vq xi^-vif / / / T« T«<^(.)T% ^(.) (a;)da^dtdsdMt;i, 

Jo Jo Jo Jo Jm-i L J 

Lebesgue's monotone convergence theorem implies 

^^^""'^Vo Jo Jo (^■)d^'d*dsdftd«i, 
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This can be written a bit shorter with potential notation 

i4i(r) -> 2Vq f x(i - yif r I u'^ \Tf^{-)z'^u'^ip(-)\ {x)dxdsdvi, 

Jo Jo Jr'^ 

Note that by assumptions (A5) the integral above is finite. 

Now we fix (5 > and choose e such that e > such that sup^g^Q — s + h) — x(l — s)| < 5 we 

[T^M-)Tns+hM-'l] i^)dxdtdsdhdvi 



Jo Jo 



Q 



It easy to deduce that the first integral is convergent (it is smaller then A^\{T) in fact). Let us deal with 
the second one. It can be written as 



F2 



(t)T« T«<^(.)T«^_^ <^(.) (x)dxdMsdM«J 
Jo Jo Jo Jo JvA '- ^ ' J g 

Let us substitute s s/T, h h/T, t t/T and recall that Ft = T^^^ 

/ / / / / l(T.,T.o('i)l(0,T{i-.o)(s)l(0,T{i-.-.O)W'7't'' p.''¥'(0^.+;.¥'(0 (a;)da;dtdsdM«i| 
./o Jo Jo Jo Jr^ m 

By assumption (A5) the integrand is dominated by integrable function 7^*^ |^7^'^^(-)7^'^^(/3(-)j (a;) hence 
Lebesgue's theorem implies the convergence to 0. We can take 5 arbitrarily small hence 

^42(7)^0. (91) 

We are left with estimation of A^iT) — Aa{T). By equation (52) and inequality (59) we have 

AiiT)- Ai{T) <2Vq f f f T,'^^[uT{-,T-s,s)vT{-,T-s,s)]{x)dsdtdx. 
Jv.-^ Jo Jo 

Using estimate (60) and (50) we write 

A4(r)-A4(r)<^ / r fTQ\rTQ^^^T{;T-s) 
Jm'' Jo Jo Uo 

Using (57), after simple calculations, we get 



(a;)dudsdtdx. 



AiiT) - AiiT) < ^ [ [ [ uT^'p{x)dudtdx. 

JvA Jo Jo 

Now, by using d'Hospital rule, it follows easily from (A6) that 

AiiT) - A^T) -* 0. 

4.2 Calculations - tightness 

We are left with proving inequalities (84) and (85). This can be done by evaluating the Ifs of the inequal- 
ities using equations derived in Section 3.4 and later estimating each of the resulting terms separately. 
Calculations are quite lengthy, for the sake of brevity in the paper we present only one illustrative example. 
Namely, consider the terms arising from the second term of (83) 

D{T) = / r r T,?^^ {{v{0)"f(x, T - w, w)) dwdhdx. (92) 
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From (80) and (81) it is easy to notice that ^"(O) < c/F^ hence 

D{T) <-^ f r f TZ^v{0)"{x, T-w, w)dwdhdx. (93) 

Jr'' Jo jo 

Now substitute v' with the first term of (81), we denote this new expression by Di (the expression 
resulting from the second term can to be estimated in a similar way) 

r-T nh 



p p 1 PIT. pW 

Di{T) = -^ / / / TZu'fT{x,T-u)vmx,T-u,u)du 

JR'' JO JO Uo 



dwdhdx. (94) 



Finally we use (81) which yields 



JWJ Jo Jo Uo \ JO J . 



Changing the order of integration and using (57) get 

^'^^^ ^ ^ Jo Jo r Jo ^^^^'''^^^^'^'''^ L^^-'^' ^^^'^ {Hx)Tu-Mx))] d^;d^.d»dMx. (96) 
Obvious changes of variables gives 



i4 fl fh fw fu 



X(l - n)x{l - V) £ T^^,_,, [t«_„, (v.(x)T«_„,^(x))] 



= / / / X{l-u)x{l-v) I T^V.-J^TV.)(vW^TV.,^W)|d«d^du,dM^. 

JO Jo JO JO 



Recall that we are using the scheme presented in Section 3.1 hence inequality (34) holds. We apply it to 
x(1 — v), use inequality T^^^_^^Lp{x) < ce"^'''"""' and integrate with respect to v 



1 /'h fW 

dudwdhdx. 



i).(T)<^(t-.)/ /'f x(i-^)£^.v.) 

Changing the order of integration and integrating with respect to w we get 

Di{T) < ^(f _ s) £ £ x(l - ^)T(h - u) £ T^^^_^^^{x)dudhdx 
Easy calculations yield 



Di{T) <c{t-s) I I x(l -u) J ^T{h- u)T^^^_^.^^{x)dudhdx. 



Using assumption (A6) one easily gets 

Di{T) < c{t-sf-'T-\ 
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